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Abstract 

By examining the structure in momentum and coordinate space of a two-body interaction spher- 
ically symmetric in its local coordinate, we demonstrate that it can be disentangled into two 
distinctive contributions. One of them is a medium-independent and momentum-conserving term, 
whereas the other is functionally -and exclusively- proportional to the radial derivative of the 
reduced matrix element. As example, this exact result was applied to the unabridged optical 
potential in momentum space, leading to an explicit separation between the medium-free and 
medium-dependent contributions. The latter does not depend on the strength of the reduced ef- 
fective interaction but only on its variations with respect to the density. The modulation of radial 
derivatives of the density enhances the effect in the surface and suppresses it in the saturated 
volume. The generality of this result may prove to be useful for the study of surface-sensitive 
phenomena. 
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I. INTRODUCTION 



During the past two decades, several developments in theoretical nuclear research have 
allowed significant improvements in the microscopic description of nuclear collisions [1, 2}. 
Such is the case of nucleon scattering, where continuous efforts have led to detailed micro- 
scopic realizations of optical model potentials, thus providing the most complete current 
first-order nonrelativistic description of nuclear collisions off nuclei over a wide energy range 
and various targets. These calculations emphasize a detailed treatment of the NN effective 
interaction, particularly its density dependence as well as its energy and momentum depen- 
dence. Quite generally, however, all realizations of the optical potential become a single 
expression in the form of a convolution of medium-dependent effective interactions folded 
with the target ground-state mixed density. In this article we demonstrate that the optical 
potential can be expressed as the sum of two very distinctive terms, one of them depends 
exclusively on the free-space t matrix and the other as a gradient of the reduced g matrix. 
This result implies that intrinsic medium effects in optical potentials become enhanced in 
the nuclear surface and suppressed in the saturated volume. 

Since the early realizations of microscopic optical potentials [3-5], the role of nuclear 
medium effects has been a major issue in the study of the dynamics involved in nucleon- 
nucleus collisions. In these studies density-dependent nucleon-nucleon local (NN) effective 
interactions have been developed to represent the force between nucleons in the nuclear 
medium. The use of these local forces, with suitable local density prescriptions, have led 
to folding optical potentials in coordinate space which provide reasonable descriptions of 
NA scattering data at energies between a few tens of MeV up to near 400 MeV. Recent 
developments [2, 6], within the same philosophy, have succeeded in including non localities 
in the optical potential stemming from the exact inclusion of the exchange term. In this case 
the full mixed density from nuclear shell models are used and provide reasonable account of 
the existing NA scattering data. 

A slightly different strategy has been followed by Arellano, Brieva and Love (ABL), with 
the realization of folding optical potentials in momentum space [7, 8]. In their approach 
genuine nonlocal g matrices, based on the Brueckner-Bethe-Goldstone infinite nuclear matter 
model, are folded to the ground-state local density of the target. As a result, nonlocal 
potentials are obtained with varying degree of success in describing the low and intermediate 
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energy data. 

The inclusion of medium effects has also been addressed within the spectator expansion 
[9], where the coupling between struck nucleons and target spectators is taken into account. 
This approach is an extension of the Watson, and Kerman, McManus and Thaler theories, 
with focus is on the many-body propagator involved in the (A+l)-body problem. 

From a more general prospective, various formal expressions of the optical potential 
can be found in the literature [10-14]. Although they may differ in the way they establish 
contact with the bare NN potential, they all become a folding expression between the target 
ground state and a generalized two-body effective interaction. In this article we analyze this 
general expression and demonstrate that, regardless of the model of utilized to represent the 
effective interaction, the intrinsic medium effects become manifest in the nuclear surface. 
The implications of this result are examined in the framework of an infinite nuclear matter 
model for the NN effective interaction. 

This article is organized as follows. In Section II we outline the general framework, 
discuss the structure of two-body operators and introduce the 'asymptotic separation' for 
spherically symmetric systems. The result is then applied to the unabridged optical potential 
in momentum space. In Section III we make use of an infinite nuclear matter model for the 
effective NN interaction and examine its implications in the optical potential. Furthermore, 
we analyze its consistency with the ABL approach, and assess the medium sensitivity of 
selected matrix elements at various energies. In Section IV we present a summary and the 
main conclusions of this work. Additionally, we have added three appendices where we 
include some intermediate steps. 

II. THE OPTICAL POTENTIAL 

A general representation of the optical model potential for collisions of a hadronic probe 
with kinetic energy E off a composite target is given by the expression 



where T represents a two-body effective interaction containing, in general, information about 
the discrete spectrum of the many-body system. The one-body mixed density p(p', p) rep- 
resents the ground-state structure of the target. Thus, a fully consistent evaluation of the 
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optical potential by means of the full f matrix would require the solution of the (A+ ^-body- 
system, a formidable task. This difficulty is circumvented by treating separately the ground 
state and the two-body effective interaction. This separation becomes suitable at interme- 
diate and high energies, where the discrete spectrum of the many-body Green's function is 
distant from the projectile energy in the continuum. Then, the target ground-state may be 
described resorting to alternative framework such as Hartree-Fock-Bogoliubov. The effective 
interaction, in turn, can be modeled using the Brueckner-Bethe-Goldstone approach. 

A. Two-body effective interaction 

Let us first focus our discussion on the two-body effective interaction and examine its 
structure. Quite generally, the representation of the two-body operator T in either momen- 
tum or coordinate space requires the specification of four vectors. We denote the coordinate 
representation of T in the form 

(rV \f\rs)= T(r's';r s) , 

where the 'prior' coordinates of each particle are r and s, respectively. Similarly, r' and s' 
refer to the 'post' coordinates of the same particles, as shown in Fig. (la). An alternative 
set of coordinates is summarized by the transformation 

R> = (r 1 + s')/2 ; x' = r' - s' ; 

(2) 

R — (r + s)/2 ; x = r — s ; 

where x represents the prior relative coordinate of the pair and R their respective center 
of mass, as illustrated in Fig. (lb). With this transformation we express the equivalence 
T(r's'; rs) = T R i R (x', x) . Thus, following the procedure outlined in Appendix A we express 
the momentum space elements T = (k'p' \ T \ k p), in the form 

f = J wr elZiw '~ w) 9Z%{W ' + w)] b ' b] ' (3) 

where gz represents the reduced interaction at the local coordinate Z. Here we have denoted 

W = k + p, b=\(k-p), 
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the prior total and relative momenta, respectively. The same construction applies to the 
post momenta, where prime marks are used. The relationship between these momenta is 
illustrated in Fig. (2). Eq. (3) for T also expresses the role of vector Z, 

Z = l(r' + r + s' + s) . 

the center of gravity of the four coordinates of the two particles. We name this the local 
coordinate, the locus where the reduced interaction is evaluated. Notice that this coordinate 
is invariant under the permutation of coordinates r ^ s. 

The above representation of the T matrix displays very clearly its dependence in terms 
of the total (W, W) and relative (6, b') momenta. Additionally, the Wigner transform in 
the R and R coordinates restricts further the structure of T, suggesting the definitions 

W ± = W'-W = (p'-p)-(k- k') ; (4a) 
W + W 

W\\ = ^— = l(k> + k+p> + p). (4b) 

The vector Wj_ represents the total momentum gained by the pair upon interaction (W = 
W + W±), whereas W\\ is the average of the prior and post total momenta. These momenta 
become perpendicular only if W' 2 = W 2 . 

B. Asymptotic separation of T 

We now examine the structure of T in the context of a finite nucleus with spherical 
symmetry. By that we understand that gz depends only on the magnitude of the local 
coordinate, \Z\ = Z. Additionally, let us assume that as Z — > oo, gz tends to its free- 
space form #00. If we omit the three vector arguments of g z and decompose it as gz = 
{gz- goo) + goo, then 

f = 5(W ± ) 9oo + -J-^JdZ e* z - w - (g z - 9o o) ■ 

Carrying out the solid angle integration, the integrand of the remaining radial integral is 
simply 4nZ 2 jo(ZW±)(gz — goo), which integrated by parts yields the asymptotic separation 

f = 5(W L ) 9oo -T^fz 3 dZ ^(Z W ± ) || . (5) 

Here S denotes the Dirac delta function and $i represents the profile function defined by 

$i(*)=ji(t)A, 
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with ji the spherical Bessel function of order 1. In Fig. (3) we plot $1, where we observe 
that its peak value (1/3) occurs at the origin, and that it is mainly contained within the 
region t < 4. 

What is interesting about Eq. (5) is that it separates unambiguously the free-space 
contribution of the T matrix from its medium-dependent counterpart. On the one hand, 
the medium dependence enters solely as the gradient of the reduced element while the total 
momentum is not conserved. On the other hand, the medium-independent contribution 
does conserve momentum, as dictated by S(W±). This contrast is physically consistent 
with our notion about non-translational invariant systems. By introducing a Z-dependent 
reduced interaction, the two-body T matrix does not expresses conservation of the total 
momentum. The conservation becomes manifest only if dgz/dZ = 0, as in the cases of 
interacting nucleons in infinite nuclear matter or free space. More interestingly, the result 
displayed by Eq. (5) is sufficiently general to allow us to model the medium dependence 
in a finite nucleus and justifies the incorporation of medium effects in distorted wave Born 
approximations (DWBA). 

The application of the asymptotic separation in Eq. (1) for U yields some undisclosed 
features. Let us first change variables from p,p' to P, Q, 

P=(p' + p)/2; Q = p'-p, 

so that dp' dp = dQ dP. These two vectors represent the mean and transferred struck- 
nucleon momenta, and the integration on them accounts for the Fermi motion of the target 
nucleons. Analogously, let us denote 

K = (k + fc')/2 ; q = k-k' ■ 

so that Wj_ = Q — q. With this notation we re-express the vector arguments of the reduced 
g matrix, 

g z {K\ h b', b) -> g z (K + P, k+) , 

where 

K ± = \[K-P±\{q + Q)]. (6) 

With these considerations and using Eq. (5) for the two-body interaction, the unabridged 
[22] optical potential takes the form 

U = U + U 1 , (7) 
6 



with 



U = J dPp(q;P) 9oo , (8a) 



t/i = - ^2 / <IQ<1P/>(Q--P) x 



^dz^zig-Afi)^. (8b) 







The first term, U , depends exclusively on the medium-free reduced matrix, whereas the 
second depends on the gradient of g. 

III. NUCLEAR MATTER MODEL 

In the preceding analysis we have made no mention to a specific approach to model the 
T matrix. In this regard, the framework is general enough to include various strategies to 
describe an effective two-body interaction in the realm of a finite nucleus. However, if the 
reduced matrix g is taken [7] as the antisymmetrized Brueckner-Bethe-Goldstone reaction 
matrix of starting energy E, 

g(E) =v + v ■ Q r g{E) , 

h + irj — hi — ti2 

then the reduced matrix at infinity, g^, corresponds to the free scattering matrix t(E). In 
the above equation hi and hi correspond to quasi-particle energies and Q, the Pauli blocking 
operator. Therefore, the first term of the optical potential in Eq. (7) becomes 

U (E) = J dP P (q;P)t(E), 

the lowest-order free t matrix full-folding optical potential in the Watson and Kerman- 
McManus-Thaler approach [12, 15]. Actual calculations of this contribution were realized 
in the early nineties [16-18]. 

With regard to U\ we note that the g matrix in Brueckner-Bethe-Goldstone approach 
is a functional of the density, g = g[p\. If the reduced matrix is evaluated at a density p 
specified by the local coordinate Z, then 

dg z _fdp\ x dg (g) 



dZ \dZ J dp 

Considering that p'(Z) = dp/dZ peaks in the nuclear surface, the intrinsic medium- 
dependent contributions to the optical potential become accentuated in that region. The 



strength of such contributions will depend on dg(E)/dp, an energy-dependent operator in 
spin-isospin space. 

To focus these ideas, in Fig. (4) we characterize the proton and neutron densities in 208 Pb 
[19], where in the upper frame we plot the proton and neutron density, in the middle frame 
the local Fermi momentum given by 

k F = (3tt 2 P ) 1/3 , 

and in the lower frame the negative gradient of the density. We have multiplied this function 
by Z 3 to account for its actual weight in the radial integration. 

What becomes clear from this figure is that medium effects stemming from dg/dp become 
dominant (if non zero) in the region between 5.5 fm and 9 fm. In this range the local Fermi 
momentum kp varies between 0.3 fm -1 and 1.2 fm -1 , suggesting the densities at which the 
g matrix needs to evaluated. 

A. Contact with the ABL approach 

The ABL approach to the optical potential constitutes an extension of the early full- 
folding approach based on the free t matrix [7]. This extension makes use of an infinite 
nuclear-matter model to represent the NN effective interaction between the projectile and 
the target nucleons. Upon the use of a simplifying assumption regarding its relative momenta 
dependence, and resorting to the Slater approximation of the mixed density, the optical 
potential takes the form of a folding of the diagonal (local) density with a nonlocal Fermi- 
averaged g matrix. We verify this result as a limit case of Eq. (8b) for the unabridged 
optical potential. 

If we assume that the g matrix in Eq. (8b) exhibits a weak dependence on the transferred 
momentum Q, then the relative momenta in g can evaluated at Q = q, consistent with the 
peak of $i at the origin. With these considerations we set 

K ± ^K^^l(K-P±q)), (10) 

and symbolize 

g z {E)^gf{E) 
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Additionally, if we use (Bl) for the Slater approximation to p(Q; P), then 

Ux = - - Z 3 dZ / Z' 2 dZ' p(Z') x 
Jo Jo 



1 



dQ ^(Z\Q - q\) j (QZ') x 

dP^^S F (P;Z'). (11) 



Here the rightmost integral does not depend on Q, therefore the dQ integration involving 
$1 j (QZ') can be performed separately. Using Eq. (C2) in Appendix C and reordering the 
integrals we obtain 

U x = -^J q Z' 2 dZ' Jo (qZ')p(Z') J dP S F (P;Z') x 

x 'rfze(z-z') dg li E) ■ (12) 



o dZ 



The integration over Z is immediate. If we identify g™(E) = t(-E), the free t matrix, then 



U x = -An l Z' 2 dZ' 3o (qZ f )p(Z f ) x 



o 



dPS F (P;Z')[t(E)-gP(E)}. (13) 

The integral involving the free t matrix leads to Uq(E), whereas the one involving gf)(E) 
leads to the ABL in-medium folding potential Uabl{E), U\ — > — Uq + Uabl, with 

POO 

U ABL {E) = An / Z' 2 dZ' jo(qZ') p(Z') gf{E) . (14) 

When this term is substituted in Eq. (7) we obtain U(E) — > Uabl(E), the expected limit. 
As seen, this result illustrates the fact that Uabl can be rigorously derived from Eq. (8b) 
using a few simplifying assumptions, i.e. Slater approximation of the mixed density and a 
weak dependence of the g matrix in the Q momentum. 



B. Contrast with the p(d/dp) rearrangement term 

During the course of this work it was noticed by some colleagues certain resemblance 
between Eq. (8b) for the unabridged optical potential, and medium corrections in the 
form of a p(d/dp) term proposed by Cheon and collaborators [20]. It becomes appropriate, 
therefore, to point out the differences in context and form of this apparent resemblance. 
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The p{d/dp) term of Cheon et al. emerges after a perturbative treatment of the transition 
density for inelastic scattering. In this case the optical potential for elastic scattering is 
given schematically by U opt = G(p)p, with G(p) the Brueckner G matrix and p the ground 
state mixed density. The transition potential for inelastic scattering, U tr , is expressed as 
U tr = ( y dU opt /dp)p tr , with p tr the transition density. Combining these equations it is shown 
that the transition potential can be expressed as 

dG( P y 



tr 



G(p)+p- 



Ptr 



dp 

As observed, the optical potential U tr results the sum of the elastic term with a corrective 
term of the form pip j dp). This correction accounts for the rearrangement of the target 
nucleons in an inelastic process. Additionally, the pip j dp) form of the correction implies 
relatively uniform contributions in the nuclear interior, in contrast with the p'p/dp) term 
in Eq. (8b) for elastic scattering, where the intrinsic medium effects manifest dominantly 
in the nuclear surface. The unabridged optical potential discussed here represents elastic 
processes and its extension to inelastic scattering would require further analysis. 



C. Medium sensitivity 

The actual evaluation of the unabridged optical potential [c.f. Eq. (8b)] constitutes a very 
challenging task beyond the scope of this work. Indeed, each matrix element requires the 
realization of a 7-dimensional integration, three more dimensions than current calculations 
in the ABL approach. However, it is possible to assess the relative importance of selected 
terms in the medium-dependent U\ contribution. In order to isolate the role of medium 
effects, let us define the amplitude 

T(Z'; E,Z) = J dP S F (P; Z') . (15) 

In the context of Eq. (12) this amplitude accounts for the Fermi average of the gradient of 
the effective interaction in the limit k± — > Thus, it is reasonable to expect that this 
quantity accounts for the leading contributions stemming from the dQ integral in Eq. (8b). 

In addition to Z', Z and E, the T amplitude depends on the momenta k and k', and spin- 
isospin degrees of freedom. To examine its radial behavior, we have evaluated the diagonal 
(Z = Z') elements which are plotted in Fig. (5) for the central pp (upper frames) and pn 
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(lower frames) channels. The real and imaginary components are shown in the left and right 
panels, respectively. These amplitudes are evaluated on-shell (E = k 2 /2m = k' 2 /2m), at 
forward angles, for nucleon energies of 65 MeV (solid curves), 100 MeV (long-dashed curves), 
200 MeV (short-dashed curves) and 300 MeV (dotted curves). The density considered in 
this analysis is that of 208 Pb, as in Fig. (4). 

By simple inspection we observe that medium effects accounted for by U\ add attraction 
and absorption to the medium-free Uq contribution. All of them are localized in the nuclear 
surface, becoming weaker in the interior. Additionally, the various components of Y exhibit 
distinctive features regarding their energy and medium sensitivity. The energy dependence 
is identified by the separation of among all four curves, being this most pronounced in the 
case of ReT pn (E) and followed by ReT pp (E). The weakest energy dependence occurs for 
the absorptive component in the pn channel. 

Regarding the intrinsic medium effects, the strongest dependence occurs for the real pn 
amplitude at 65 MeV. In contrast, the weakest dependence occurs for the absorptive pp 
(upper-right frame) amplitude. Additionally, with the exemption of ImT pn (E), all the other 
amplitudes exhibit decreasing strength with increasing energy. 

IV. SUMMARY AND CONCLUSIONS 

We have examined the role of medium effects in the optical model potential for hadron- 
nucleus scattering. The analysis is based on a close scrutiny of the structure of the two- 
body effective interaction with spherical symmetry in its local coordinate Z, leading to its 
asymptotic separation. As a result, we demonstrate that the unabridged optical potential 
can be separated into two very distinctive contributions. One of them is a momentum- 
conserving and medium-independent term, while the other is functionally proportional to the 
radial derivative of the reduced matrix element. If the Brueckner-Hartree-Fock g matrix is 
used to model the NN effective interaction, then the medium-independent term of the optical 
potential corresponds to the well known Watson-KMT lowest order full-folding optical model 
potential, while the medium-dependent term depends exclusively on the gradient of the 
reduced g matrix. The modulation by the radial derivatives of the density enhances that 
effect in the nuclear surface and suppresses it in the saturated volume. 

The assessment of the intrinsic medium effects by means of the T(Z';E, Z) amplitude 
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points to stronger medium sensitivity in the real part of the pn amplitude at 65 MeV, in 
contrast with the absorptive pp component. The energy-dependence of these effects are 
relatively weak in the absorptive pp and pn channels. The introduction of the T amplitude 
may prove to be very useful for the assessment of intrinsic medium effects, such as those 
associated with the inclusion of non spherical components of the Pauli blocking studied in 



Although we have not provided a full realization of the unabridged optical potential, we 
have been able to extract its exact functional dependence in terms of the nuclear density 
when the system is spherically symmetric. This feature may become useful in the context 
of semi-phenomenological approaches, where the free-space contribution may well be repre- 
sented in terms of impulse-approximation-like potentials, while the medium dependent term 
can be modeled as function of p'd/dp couplings. At a more fundamental level, it would be 
interesting to identify missing features in the ABL approach relative to the unabridged po- 
tential. As inferred from Fig. (5), intrinsic medium effects become most pronounced -in the 
nuclear surface- at lower energies. This surface- sensitive phenomenon may be of particular 
importance in the study of rare isotope beams, where highly unstable nuclei are collided 
against hydrogen targets. In this case, the traditional intermediate energy regime is reached 
with 30A-100A MeV beams. 

APPENDIX A: REPRESENTATION OF QUANTUM TWO-BODY OPERA- 
TORS 

Let A a two-body operator and denote its coordinate representation by (rV | A \ r s) — 
A(r's'; r s) , where r, s (r', s') represent the prior (post) coordinates of each particle. An 
alternative set of coordinates is summarized by Eqs. (2) and illustrated in Fig. (lb). With 
this transformation we can express A(r's'; rs) = A R/R (x', x) . If we denote A = {k'p' | A | 
k p) , then 



Ref. [21]. 
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In terms of the coordinate set defined in Eqs. (2) we re-express the above integral in the 
form 



A = - — — / dR' dR dx' dx x 

(2tt) 6 J 

e -i { W<.B! + V. x >-W.R- b . x ) ARiR{x ^ x) ) 

where we denote W = k + p, the momentum of the pair prior interaction and b = (k — p)/2 : 
the relative momentum. Here again prime marks denote post interaction as shown in Fig. 
(2). To the above expression we apply the Wigner transform in the prior and post center of 
mass coordinates R and R'. Defining Y = R' — R, and Z = (R' + i?)/2 we obtain 

A = J WT etZ<W '~ W) 9z[\{W + W); b>, b] , (Al) 
where gz represents a reduced interaction g at the coordinate Z and is given by 
g z (K li; b',b) = J ^e- 1 ™ 11 J dx'dxe-W^'-^x 



APPENDIX B: THE GROUND-STATE MIXED DENSITY 



Let us denote the one-body density matrix in momentum space by p(p f , p). In coordinate 
space this matrix is represented by p{r',r), where 

p(p',p) = j^y 3 J dr'dr e"* ' r ' p(r< , r) e^ r . 

It is customary to introduce the center-of-mass Z = (r' + r*)/2, and relative x = r' — r, 
coordinates. Similarly, we define the mean P = (p f + p)/2, and transferred Q = p' — p, 
momenta. With the use of these transformations we denote p(p',p) = p(Q; P) ', p(r', r) = 
p(Z;x) . Therefore, 

P(Q; P) = j^JdZdr e^ p - x+ ^p{Z- x) ■ 

" ^fS dZe ' lQ ' zp{z)G{Z]Ph 



(Bl) 



where we have defined 



p{Z) G(Z; P) = J dx e- iP x p(Z; x) 
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In the Slater approximation for the mixed density its coordinate representation takes the 
form 

p(Z;x) = p(Z)F(Z;x), (B2) 

with F(Z\ x) — 3 ji(kzr) / (kzr) . Here kz, the local Fermi momentum, depends on the local 
density p(Z) through 

k = (3tt 2 P ) 1/3 • 

Within this approximation, and assuming spherical symmetry in the local density, the func- 
tion G(Z; P) can be calculated directly. When used to evaluate Eq. (Bl) for p(Q; P) we 
obtain 

POD 

p(Q; P) = 4vr / Z 2 dZ 3o (QZ) p(Z) S F (P; Z) , (B3) 
Jo 



where 



S F (P; Z) 



\^k\ 



Q{k z - P) 



(B4) 



APPENDIX C: EVALUATION OF j dQ j $i 



We evaluate the volume integral £l q (Z', Z) defined by 

where we express \Q — q\ = \/Q 2 + q 2 — 2Qqu, with u — Q-q. The solid-angle integration 
is straightforward, leading for Q 



n q (z',z) = 



4ir 
Z 2 Z' 



Q sinQZ'dQ 
q 2 -Q 2 



x 



cos(QZ) sin(gZ) sin(<5^) cos(gZ) 



qZ QZ 

This integral can be evaluated analytically and the expression depends on the location of Z 
relative to Z' . The two cases of interest are Z > Z', and Z < Z', where we obtain 

if Z < Z'\ 

j Q {qZ') if Z > Z>- 

When Z — Z' the result is finite but discontinuous. Due to its involved form we have 
preferred to omit it for clarity, keeping in mind that at Z — Z' its contribution as integrand 
vanishes. Thus, without loss of generality we can express 



n q (z',z) 



(Cl) 



2tt 2 

Q q (Z',Z) = — j (qZ') xQ(Z-Z'), 



(C2) 
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with Q(x) the Heaviside step function. 
ACKNOWLEDGMENTS 

H.F.A. acknowledges partial support provided by FONDECYT under grant 1040938. 



[1] L. Ray, G. W. Hoffmann and W. R. Coker, Phys. Rep. 212, 223(1992). 

[2] K. Amos, P. J. Dortmans, H. V. von Geramb, S. Karataglidis, and J. Raynal, Adv. in Nucl. 

Phys. 25, 275 (2000). 
[3] J. P. Jeukenne, A. Lejeune and C. Mahaux, Phys. Rep. 25, 83(1976). 
[4] F. A. Brieva and J. R. Rook, Nucl. Phys. A291, 317 (1977); ibid. A307, 493 (1978). 
[5] H. V. von Geramb, in The Interaction Between Medium Energy Nucleons in Nuclei, edited by 

H. O. Meyer (AIP, New York, 1983); 
[6] J. Raynal, computer code DWBA98 (NEA 1209/05, 1999). 
[7] H. F. Arellano, F. A. Brieva, and W. G. Love, Phys. Rev. C 52, 301 (1995). 
[8] H. F. Arellano and H. V. von Geramb, Phys. Rev. C 66, 024602 (2002). 
[9] C. R. Chinn, Ch. Elster, R. M. Thaler, and S. P. Weppner, Phys. Rev. C 52, 1992 (1995). 
[10] K. M. Watson, Phys. Rev. 89, 575 (1953). 

[11] H. Feshbach, Ann. Phys. (N.Y.) 5, 357 (1958); ibid. 19, 287 (1962). 

[12] A. K. Kerman, H. McManus, and R. M. Thaler, Ann. Phys. 8, 551 (1959). 

[13] A. L. Fetter and K. M. Watson, in Advances in Theoretical Physics, Vol. 1, edited by KA. 

Brueckner, (Academic Press, New York, 1965) 
[14] F. Villars, in Fundamentals in Nuclear Theory, edited by A. de-Shalit and C. Villi, (IAEA, 

Vienna, 1967). 

[15] H. F. Arellano, F. A. Brieva, and W. G. Love, Phys. Rev. Lett. 63, 605 (1989). 
[16] H. F. Arellano, F. A. Brieva, and W. G. Love, Phys. Rev. C 41, 2188-2201 (1990) 
[17] Ch. Elster, T. Cheon, E. F. Redish and P. C. Tandy, Phys. Rev. C41, 814 (1990). 
[18] R. Crespo, R.C. Johnson and J.A. Tostevin, Phys. Rev. C41, 2257 (1990). 
[19] J. -F. Berger, M. Girod, and G. Gogny, Comput. Phys. Commun. 63, 365 (1991). 



15 



[20] T. Cheon, K. Takayanagi, and K. Yazaki, Nucl. Phys. A437, 301 (1984); ibid. A445, 227 
(1985). 

[21] E. J. Stephenson, R. C. Johnson, and F. Sammarruca, Phys. Rev. C 71, 014612 (2005) 
[22] By unabridged we mean the nine-dimensional integration comprised of six dimensions in mo- 
mentum space (dPdQ) and three in coordinate space (dZ). Spherical mass distribution re- 
duces the dimension of non trivial integrals to seven (6+1). 







1 

x' 


i 

X 




x' 


Y 




s' 






R 






XT 



(a) ' (b) (c) 



FIG. 1: Schematic representation of the vector coordinates in a two-body operator. 




FIG. 2: Schematic representation of the momenta in a two-body operator. 
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FIG. 3: The profile function &i(t) as function of t. 
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FIG. 4: Local density p (upper frame), local momentum kp (middle frame) and Z 3 times the 
negative gradient of the density (lower frame) as functions of the local coordinate Z. The solid 
and dashed ci 




1 1 


1 1 1 1 1 

Im 


' — 1 — — h 






— | 1 1 1 

Im 














, 1 


i,i,i 



10 2 4 6 

Z [fm] 



8 10 



FIG. 5: Fermi integral of the gradient of the g matrix (Tnn) as function of the radial distance Z 
in 208 Pb. The pp (pn) channel is shown in the upper (lower) frame, whereas the real (imaginary) 
components are plotted in the left (right) panels. The solid, long-dashed, short-dashed and dotted 
curves represent results for Tnn at nucleon energies of 65, 100, 200 and 300 MeV, respectively. 
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